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Abstract 

We present a theory for the kinetics of surfactant adsorption at the interface 
between an aqueous solution and another fluid (air, oil) phase. The model relies 
on a free-energy formulation. It describes both the diffusive transport of surfactant 
molecules from the bulk solution to the interface, and the kinetics taking place at the 
interface itself. When applied to non- ionic surfactant systems, the theory recovers 
results of previous models, justify their assumptions and predicts a diffusion-limited 
adsorption, in accord with experiments. For salt-free ionic surfactant solutions, elec- 
trostatic interactions are shown to drastically affect the kinetics. The adsorption 
in this case is predicted to be kinetically limited, and the theory accounts for un- 
usual experimental results obtained recently for the dynamic surface tension of such 
systems. Addition of salt to an ionic surfactant solution leads to screening of the 
electrostatic interactions and to a diffusion-limited adsorption. In addition, the 
free-energy formulation offers a general method for relating the dynamic surface 
tension to surface coverage. Unlike previous models, it does not rely on equilibrium 
relations which are shown in some cases to be invalid out of equilibrium. 
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1 Introduction 



Aqueous solutions of surface-active agents (surfactants) play a major role in various fields 
and applications, such as biological membranes, petrochemical processes, detergents, etc. 

In some important cases, equilibrium properties of the surfactant adsorption at in- 
terfaces are not sufficient, and knowledge of the kinetics is required. Processes of fast 
wetting, foaming and stability of thin soap films may serve as good examples. The kinet- 
ics of surfactant adsorption have been addressed by experimental and theoretical studies 
since the 1940s, and various experimental techniques have been devised, primarily aimed 
at the measurement of dynamic interfacial tensions 0. 

The pioneering theoretical work of Ward and Tordai ||^ considered a diffusive transport 
of surfactant molecules from a bulk surfactant solution to an interface and formulated the 
time-dependent relation, 



2cbVi- C -^^^^=dT 



/o \/t-T 

where Cb is the bulk concentration and D the surfactant diffusivity. This equation gives 
only one relation between the surface density of surfactants adsorbed at the interface, 
and Ci(t), the surfactant concentration at the sub-surface layer of solution. 

Subsequent theoretical research has focused on providing the second closure relation 
between these two variables by introducing a certain adsorption mechanism at the in- 
terface. Various relations have been suggested, resembling equihbrium isotherms 0-0, 
or having a kinetic differential form [0-[|lO|- Such theories have been quite successful in 
describing the experimentally observed adsorption of common non-ionic surfactants. Yet, 
they suffer from several drawbacks: (i) the closure relation between the surface density 
and sub-surface concentration, which expresses the kinetics taking place at the interface, 
is introduced as an external boundary condition, and does not uniquely arise from the 
model itself; (ii) the calculated dynamic surface tension relies on an equilibrium equation 
of state, and assumes that it also holds out of equilibrium [|ll|]; (iii) similar theories can- 
not be easily extended to describe more complicated systems, such as ionic surfactant 
solutions llT2i. 
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In the current work we would like to present an alternative approach to the kinetics of 
surfactant adsorption, overcoming these drawbacks. In Sec. |^ we lay the foundations of 



our model |13[, based on a free-energy formulation. Non-ionic surfactants are considered, 
recovering results of previous models and justifying their assumptions. In particular, we 
show that the adsorption of non-ionic surfactants is in general limited by diffusion from 
the bulk solution. In Sec. ^ we modify the theoretical framework and apply it to salt-free 
ionic surfactant solutions. A few models have been proposed for describing the kinetics 
of ionic surfactant adsorption ||T^-|[1^, yet none of them is able to account for recent 
experimental results for the dynamic surface tension of salt-free ionic surfactant solutions 



. We show that the adsorption in such systems is limited by the kinetic processes at the 



interface. Consequently, we point out a problem common to all previous models. Using 
our model, we then account for the recent experimental findings. Section ^ considers ionic 
surfactant solutions with added salt. The adsorption is shown in this case to be limited 
again by diffusion, and the effect of salt concentration is examined. Finally, we present a 
few concluding remarks in Sec. || and point out possible future prospects. 

2 Non-Ionic Surfactants 

Consider an interface between an aqueous solution of non-ionic surfactants and an air or 
oil phase. The system is schematically illustrated in Fig. 1. We assume that the width 
of the interface is much smaller, and its radius of curvature much larger, than any length 
scale relevant to the adsorption process. Hence, the interface can be regarded as sharp 
and flat, lying at the plane x = 0, and the problem is reduced to one dimension. At 
X — s> oo, the solution is in contact with a bulk reservoir of surfactant molecules, where the 
chemical potential and surfactant volume fraction are fixed to be fib and (pb, respectively. 
We consider a dilute solution, i.e., the surfactant volume fraction is much smaller than 
unity throughout the solution. The concentration is also smaller than the critical micelle 
concentration (cmc), so the surfactants are dissolved only as monomers. At the interface 
itself, however, the volume fraction may become large. 
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2.1 The Free Energy 

We write the excess in free energy per unit area due to the interface [i.e., the change 
in interfacial tension), A7, as a functional of the surfactant volume fraction in the bulk 
solution, > 0), and its value at the interface, (po, 

POO 

A7[0]= / A/[0(x)]rfx + /o(0o). (2.1) 
Jo 

The first term is the contribution from the bulk solution. A/ being the excess in free 
energy per unit volume over the bulk, uniform state. The second is the contribution 
from the interface itself, where /o is the free energy per unit area of the surfactant at the 
interface. The sharp, "step-like" profile considered, has led us to treat the bulk solution 
and the interface as two coupled sub-systems, rather than a single one [l^ . 

The bulk sub-system is considered as an ideal, dilute solution, including only the ideal 
entropy of mixing and the contact with the reservoir and neglecting gradient terms, 

A/(0) = l{T[01n0 - - (0,ln0, - 0,)] - /x,(0 - 0,)}, (2.2) 
a"* 

where a denotes the surfactant molecular dimension and T the temperature (we set the 
Boltzmann constant to unity). 

At the interface, however, since 0o may become much larger than (j){x > 0), we 
must take into account the finite molecular size and the interactions between surfactant 
molecules, 

M(t>o) = ^ |T[0oln0o + (1 - 0o) ln(l - 0o)] - «0o - - /^i0o| • (2.3) 

The term in the square brackets is the entropy of mixing, this time in its complete form, 
since 0o is not necessarily small. The second accounts for the energetic preference of 
the surfactants to lie at the interface, a being positive by the definition of our molecules 
as surface-active. The third is the energy of lateral interaction between neighboring 
surfactants at the interface, where /? is assumed to be positive too, i.e., expressing an 
overall attractive interaction. The last term accounts for the contact with the solution 
adjacent to the interface, where the chemical potential is fii = fi{x —>■ 0) 
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Variation of A7 with respect to 0(x) yields the excess in chemical potential at a 
distance x from the interface, 



A/i(x) 
A/io 



2 5A7 



Tln0(x) — /ife ; X > 



Tin- 



- a - /:((/)o - /ii. 



(2.4) 
(2.5) 



From Eq. ( p.4| ) we can deduce, as expected, 



/ife = T In 0b 
/ii = Tin 01, 



(2.6) 



where (pi = (j){x 0) denotes the surfactant volume fraction at the sub-surface layer. 



2.2 Thermodynamic Equilibrium 



In equihbrium, the chemical potential is equal to fib throughout the entire system (the 
variations of A7 vanish). From Eq. (p.4| ) we obtain the equilibrium profile. 



z > 0, 



(2.7) 



and from Eq. 



the equilibrium adsorption isotherm. 



(2.^ 



We have recovered the Frumkin adsorption isotherm, which reduces to the well known 



Langmuir adsorption isotherm [|l9l when the interaction term is neglected {(3 = 0). From 
Eqs. (|2.3| ), ( p.5|) and (|2.7| ) one also obtains the equilibrium equation of state, 



A7 



1 



Tln(l 



^0) + ^<Pl 



(2.9) 



which was previously derived from other, though equivalent considerations (integration 
of the Gibbs equation) B . 
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2.3 Out of Equilibrium 



Throughout our analysis we assume proportionahty between velocities and the potential 
gradient pO[, and take the surfactant mobility to be D/T according to the Einstein 
relation {D being the surfactant diff usivity) . At positions not adjacent to the interface 
this leads to the following surfactant current density, 



■( \ ^^^^ n^'^ 



(2.10) 



Applying the continuity condition, d(j)/dt = —dj/dx, we get the ordinary diffusion equa- 
tion, 

rifh rP'rh 

(2.11) 



d(j) _ jjd^(l> 
dt dx'^ 



The proximity of the interface requires a more careful treatment. First, we discretize 
expression ( |2.1| ) on a lattice with cells of size a, 

oo 

A7[0] = a5:A/(0O + /o(0o), 

i=l 

where (pi = (j){x = ia). Discretized current densities, ji, can be similarly defined. Since 
we do not allow molecules to leave the interface towards the other (air, oil) phase (i.e., 
jo = 0), we have from the continuity condition 

d(j)o ji 



and can therefore write 



dt 



dt 



32-31 



a 



D dcj) 
a dx 



d(j)o 
dt 



[2.12) 



Applying the Laplace transform to Eqs. ( |2.11D and ( |2.12D while assuming an initial uniform 
state, 0(x, t = 0) = 0b, a relation is obtained between the surface and sub-surface volume 
fractions, 0o and 0i, 



a V TT 



2(j)bVt 



\Jt-T 



dr 



+ 20b - 01. 



(2.13) 



This relation is similar to the classical result of Ward and Tordai 0], Eq. ( |1 . 1| ) , except 
for the term 20^ — 0i. This difference is due to the fine details we have considered near 
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the interface and our initial condition. Ward and Tordai's analysis assumes a continu- 
ous profile up to the non-aqueous phase, and hence replaces Eq. ( p. 121) with a simpler 
condition, dcpo/dt = {D/a)d(j)/dx\x=o- In addition, it requires an initial empty interface 
[(t>o(t = 0) = 0], whereas we set 0o(^ = 0) = 0^. At any rate, the difference vanishes when 
the cell size, a, goes to zero. Finally, we find the equation governing the kinetics at the 
interface itself pT |, 



d(f)o 



a 



In 



1 - 



a pd) 
H h — 

' rj-\ ' rj-} 



(2.14) 



dt a a T a 

Simultaneous solution of Eqs. ( p.l3| ) and (|2.14|) yields the solution of the adsorption 
problem, i.e., the time-dependent surface coverage, 0o(^)- 



2.4 Limiting Cases for the Adsorption 

In writing the above equations, we have separated the kinetics of the system into two 
coupled kinetic processes. The first takes place inside the bulk solution and is described 
by Eqs. ( |2.11|) and ( p.l2| ) [or, alternatively, by Eq. ( p.l3|) ], whereas the second takes place 



at the interface and is described by Eq. ( p.l4|) . Two important limiting cases correspond 
to the relative time scales of these two processes: 

• Diffusion-limited adsorption (DLA) applies when the equilibration process inside 
the solution is much slower than the one at the interface. One can then assume 
that the interface is in equilibrium at all times with the adjacent solution, i.e., 



the variation ( |2.5| ) vanishes, and 0o immediately responds to changes in 0i via the 



equilibrium isotherm. 

• Kinetically limited adsorption (KLA) takes place when the kinetic process at the 
interface is the slower one. In this case, the solution is assumed to be at all times in 
equilibrium with the bulk reservoir, i.e., the variation (|2.4|) vanishes, and 0o changes 
with time according to Eq. ( |2.14| ). 

One may suggest an alternative way of looking at the same limiting cases, the use- 
fulness of which will become evident later on. Let us re-examine the expression for the 
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interfacial contribution, /o [Eq. ( |2.3| )]. The DLA case corresponds to the following descrip- 
tion (see Fig. 2a). The interface is all the time at the minimum of the curve /o(0o); yet 
the shape of the curve changes with time as /ii is changed by diffusion, until it attains the 
value of /ift. The surface coverage increases with time as the minimum of /o is shifted to 
larger values of 0o- On the other hand, KLA corresponds to a different scenario (Fig. 2b). 
The shape of /o(0o) remains fixed since fii is constantly equal to fib- The surface coverage, 
(po, increases until finally reaching the value corresponding to the minimum of Jo- 
in order to figure out whether one of the above limits applies to non-ionic surfactant 
adsorption, the time scales of these two limiting cases must be compared. Let us start 
with the DLA case and look for the asymptotic time dependence of the process. We return 
to the Laplace transform of Eqs. ( p.ll| ) and ( |2.12| ) and let the conjugate variable of the 



transform approach zero |^ . In addition, since the DLA limit is currently considered, we 
neglect the kinetics at the interface and take 0o to be almost constant. After inverting back 
to t-space, this procedure leads to the following asymptotic time dependence, previously 
mentioned by Hansen p. 



O'4'0,eg 



t^OO, (2.15) 



where 0o,e(j is the equilibrium surface coverage. Looking at Eq. ( p.l5|) we can identify the 
characteristic time scale of diffusion, 

Typical values of a'^/D correspond to very short times (of order 10~^ sec), but the prefactor 
of {(po^eq/(pbY is usually very large (of order, say, 10^^). Thus, the diffusion time scales 
may reach minutes. 

In order to estimate the time scale of the kinetics taking place at the interface, we 
look at the asymptotic behavior of Eq. (p.l4|) and find 



r, ~ fK^V ^^-{a+f5^o,.,)/T_ ^2.17) 
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Since the value of D at the interface is not expected to be drastically smaller than that 
inside the solution, comparison of Eqs. ( p.l6[) and (|2.17|) leads to the conclusion that 
Trf > Tfc- We thus expect, in general, that non-ionic surfactants should exhibit diffusion- 
limited adsorption. This is one of the conclusions of our study as applied to non-ionic 
surfactants, and has been observed for quite a large number of such surfactants 0, . 
It is somewhat expected, since we did not include in the interfacial free energy, Eq. ( |2.3|) , 
any potential barrier that might lead to kinetic limitations. The "footprint" of DLA is 
the asymptotic time dependence (see Eq. ( |2.15|) and Ref. 0) 

For any reasonable dependence between the surface tension and surface coverage, 7(t) = 
7[0o(^)]; we expect the surface tension to also exhibit a similar asymptotic time depen- 
dence, 7(t) — ~ y^Trf/t. It should be compared with the asymptotic exponential time 
dependence in the case of KLA [Eq. ( |2.17| )]. Four examples of experimental results are 
given in Fig. 3, all exhibiting an asymptotic t^^/^ behavior, as expected for DLA. 

Having realized that the adsorption is diffusion-limited, the mathematical solution of 
the problem amounts to the simultaneous solution of two equations: (i) the Ward-Tordai 
equation, (|2.13|) , accounting for the diffusive transport; (ii) an isotherm, such as (|2.8|) 
(with 01 replacing 0f,), describing the (immediate) response of the surface coverage to 
changes in the sub-surface layer of solution. Useful numerical schemes were suggested for 
performing such a task 0, yielding the time-dependent surface coverage, (po{t). 

2.5 Dynamic Surface Tension in a DLA Process 

Since most experiments measure dynamic surface tensions and not surface coverages, we 
still need a relation between these two variables in order to relate theoretical calculations 
to actual measurements. We return, therefore, to the evolution of the interfacial tension 
during a DLA process. As was stated above, in this limit the interfacial contribution, 
/o(0o)) is all the time at its minimum. We can, therefore, write 

A7M= r Af[<P{x)]dx + \ 
Jo 



Tln(l-0o) + |0^ 
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If, in addition, the contribution from the bulk solution is neglected (recalling that it 
completely vanishes when equilibrium is reached), we are left with the equilibrium relation, 
Eq. ( p.9|) . Hence, the equilibrium equation of state, relating the surface tension to the 
surface coverage, holds approximately also out of equilibrium. Note, that this conclusion 
is valid only in the case of diffusion-limited adsorption. 

Practically all previous works assumed that the equilibrium relation between the sur- 
face tension and surface coverage holds for the dynamic surface tension as well. As we 
have just found, non-ionic surfactants usually do undergo DLA. Hence, the assumption 
employed by previous works was justified, as far as non- ionic surfactants were concerned. 
Indeed, satisfactory agreement with experimental findings is obtained, when results of the 
numerical schemes mentioned above are related to the dynamic surface tension via the 
equilibrium equation of state However, this conclusion is drastically modified for 

ionic surfactants, as we show in the next section. 

The dependence of A7 on 00, as defined by Eq. (p.9| ), is shown in Fig. 4a. Note 
the moderate slope in the beginning of the process; the surface coverage significantly 
changes without a corresponding change in the surface tension. It is a result of the 
competition between the entropy and interaction terms in Eq. (|2.9|) . As the surface 
coverage increases, the surface tension starts falling until reaching its equilibrium value. 
Since 0o monotonically increases with time during the adsorption, we expect the time 
dependence of A7 to resemble the schematic curve depicted in Fig. 4b: a slow change in 
the beginning, then a rapid drop, and eventually a relaxation towards equilibrium. This 
is, indeed, in agreement with dynamic surface tension measurements {e.g., 0). Returning 
to Fig. 4a, the surface tension will start its rapid fall roughly when the second derivative 
of A7 with respect to 0o changes sign, i.e., when 

As one examines surfactant solutions of increasing bulk concentrations (but always below 
the cmc), the surface coverage corresponding to the beginning of the drop in surface 
tension will be reached earlier along the process. The initial period of slow change in the 
tension will shrink, until finally vanishing behind the finite experimental resolution. This 
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trend is observed experimentally 

The need for an interaction between surfactant molecules in order to account for 
such a time dependence of surface tensions was previously realized by Lin et al. . 
They even suggested the existence of a transition from a gaseous to a liquid phase as 
being responsible for the initial period of almost constant tension. In fact, the form 



of /o as defined in Eq. (2^) may account for a two-phase region, but only if /5 > AT. 
As demonstrated above in Fig. 4 (where (3 = ST), this is not necessary, however, for 
recovering the experimentally observed time dependence. 



3 Ionic Surfactants without Added Salt 

We now consider the problem of ionic surfactant adsorption. The main difference com- 
pared to the previous, non-ionic case is the introduction of electrostatic interactions. The 
kinetics of the system include, apart from the diffusive transport of molecules and their 
adsorption at the interface, the formation of an electric double layer due to the increasing 
surface charge. We start with the case of a salt-free solution, where the only charges 
present are those of the surfactant ions and their balancing counter-ions. The system 
is schematically illustrated in Fig. 5. In such a case the electrostatic interactions are 
unscreened, and thus, as we shall see, have a drastic effect on the adsorption process. 
Instead of the single degree of freedom we have specified in the non-ionic case, namely 
the surfactant profile, 0(x,t), we should consider here three degrees of freedom: the sur- 
factant ion profile, t), the counter-ion profile, (f)~{x,t), and the local mean electric 
potential, tp{x,t) (without loss of generality, we take the surfactant ions as the positive 
ones) . 

For simplicity, we assume the following: (i) the surfactant molecules are fully ionized; 
(ii) the surfactant ions and counter-ions are monovalent (extension to general valencies, 
however, is quite straightforward within our mean- field formulation); (iii) as in the previ- 
ous section, the solution is assumed to be dilute. In the current context this last assump- 
tion also allows us to treat the surfactant and counter-ion volume fractions as independent 
variables, as the probability for an ion and a counter-ion to overlap at the same spatial 



11 



position is negligible. In addition, and without loss of generality, we take the surfactant 
ions to be the positive ones. 

The extension of the model to ionic surfactants evidently demonstrates the benefits of 
the free-energy formulation we have employed. We just repeat the scheme of Sec. ^ while 
adding terms for the counter-ions and electrostatic interactions. Unlike previous models 
and due to its simplicity, this formulation will allow us to clarify the complex problem of 
ionic surfactant adsorption and reach novel conclusions concerning it. 



3.1 The Free Energy 

Following the analysis of Sec. ^, we write the change in interfacial tension in the form 

POO 

A7[0+,0-,^]= / {A/+[0+]+Ar[0-] + /,,[0+,0-,^]}cia; + /o(0+,^o). (3.1) 



The first two terms are the contributions from the bulk solution, depending on the surfac- 
tant and counter-ion volume fractions. The third term is the electrostatic energy stored 
in the bulk solution. The last term is the contribution from the interface itself, depending 
on the interfacial values of the surfactant volume fraction, 0,}, and the electric potential, 
ipQ. The counter-ions are assumed to be surface- inactive, and their contribution to this 



term is neglected As in the previous discussion, the two types of ions dissolved in 
the bulk solution have contributions coming from their entropy of mixing and chemical 
potentials, 

A/^(0^) = ^{T[0±ln0± - 0± - (0^1n0± - 0±)] - - 0±)}, (3.2) 

where are the molecular dimensions of the two ions and 0^, /i^ their volume fractions 
and chemical potentials, respectively, at the bulk reservoir. Our assumption of monovalent 
ions implies 

_ <f>b _^ 



(a+)3 (a-)3 
where Cfe is the bulk concentration. 



12 



The electrostatic term contains the interaction between the ions and the electric field 
and the energy associated with the field itself, 

where e is the electronic charge and e the dielectric constant of the solvent (water). Finally, 
the modified expression for the interfacial contribution, /q, is obtained from Eq. (|2.3| ) just 
by adding an electrostatic term, 

foi<f>o. ^o) = ^ |t[0+ In 4 + (1 - 0+) ln(l - 0+)] - - §(0o+)' - f^t<Pt + e<l>t^o 

(3.4)^ 

When the interface is still almost uncharged, another electrostatic contribution should 
be considered, namely that of the interaction between the ions and their "image" charges 
beyond the interface, as was pointed out by Onsager and Samaras . This force decays 
like and since the dielectric constant of the other phase (air, oil), is smaller than 
that of the solvent, e, it is repulsive. The "image" force is comparable to the repulsion 
from the surface charge for distances 



/ / ^ e-eo + 

X < w X Ci ■ 

V 167r0+ e + e„ 

Obviously, when the right-hand side of this expression becomes smaller than the molecular 
dimension, a"*", this force is irrelevant. For a water- air or water-oil interface this happens 
when (pQ ~ 0.02, i.e., very soon along the adsorption process. We allow ourselves, there- 
fore, to disregard such early stages of adsorption and neglect this interaction altogether 
(it should be mentioned, however, that in the case of surface- macizt;e electrolytes, the 
"image" interaction has a significant contribution to the interfacial tension). 

We now take the variation of A7 with respect to (/)^(x) to get the excess in electro- 
chemical potentials, 

A^^(x) = fx^{x) - fj.^ = {a^Y^^ = Tln^^ ±etp - fxf ; x>0 (3.5) 

^/^o = f^o - f^i = (a'^y^^ = T In a- /3(f)^ + e'^o- f^t- (3.6) 

0(1)0 1-00 
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We require that the electric potential vanish far away from the interface, that is, iplx 
oo) = 0, and hence, deduce from Eq. (|375| ), as expected. 



= Tln0± 

fif = Tln0f±e?/'i, (3.7) 

where (pf and ipi denote, respectively, the sub-surface values of the surfactant volume 
fraction and electric potential. When we take the variation of A7 with respect to tp{x) 
and set it to zero (since only electrostatic effects are considered) , the Poisson equation is 
obtained, 

^^_4vre/^+ 0_\ 

dx^ e \{a+f {a-fj' ^ ' ' 

Finally, variation with respect to ip{x = = ipQ yields the expected boundary condition. 



dil) 
dx 



~ 6(a+)2^°' 



- (3.9) 



which is equivalent to the requirement of overall charge neutrality (only surfactant ions, 
and no counter-ions, are allowed to be adsorbed at the interface). 



3.2 Thermodynamic Equilibrium 

In equilibrium, the variations ( p.5|) and (|3.6|) vanish, and we recover the Boltzmann dis- 
tributions, 

0±(x) = (/.^e^^'^^")/^; X > 0, (3.10) 



and the adsorption isotherm. 



+ Q-{a+f3,j>+-eiPo)/T 



(3.11) 



We have recovered the Davies adsorption isotherm for ionic surfactants Combining 
Eqs. ( p.8|) and ( p.lOj ) leads to the well known Poisson- Boltzmann equation, 

^ = ^smh-, (3.12) 
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which determines the equihbrium double- layer potential |3^. Integrating this equation 
once and using the boundary condition yield the surface potential, 



2T 

— sinh-^(A0+), 
e 



(3.13) 



where 



A 



and K ^ = [Sti CbC^ / eT) ^1"^ is the Debye-Hiickel screening length. Substituting i/jq from 
Eq. ( 3.13 ), the Davies isotherm can be expressed in terms of alone, 

= , (3.14) 

0+ + [A0+ + ^(A0o+)2 + l]^e-(-+f^^o)/T 

From Eqs. ( |3.10|) , ( p.l2| ) and ( ^lOf ) one can calculate the contribution of the bulk solution 
to the interfacial free energy at equilibrium. 



req 
J bulk 



/ {Af^ + Af- + U)dx 
Jo 



2T 



(a+)^X 



cosh — 1 
2T 



2T 



(a+)2A 



[v'(A0+)2 + l-l]. (3.15) 



From Eqs. (|3.4|) , (|3.6| ) and (|3.15|) we then get the equilibrium equation of state, 



A7 



Tln(l 



2T 



\/(A0o+)^ + 1 



(3.16) 



In the limit of weak electric fields {X(f)Q <^ 1), the electrostatic repulsion between sur- 
factant ions at the interface predominates. As a result, the electrostatic correction to 
the non-ionic equation of state ( |2.9[ ) is quadratic in (pQ and effectively reduces the lat- 
eral attraction term. In the limit of strong electric fields {X4>q ^1), however, the high 
concentration of counter-ions near the interface makes the electrostatic term become only 
linear in (ht . 



3.3 Out of Equilibrium 

We write the current densities as in the previous section. 
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where are the diffusivities of the two ions. Applying the continuity condition, the 
Smoluchowski diffusion equations are obtained. 



d 



dt 



dx \ dx 



dx, 



(3.18) 



As in Sec. ||, we treat separately the positions adjacent to the interface by discretizing 
the expressions for the various contributions and considering the current densities near 
the interface. The condition (see Sec. p73|) 

Ji+ 



dt 



leads in this case to 



D 



± 



dt 



dx 



dip 



T^^ dx 



dt 



(3.19) 



The kinetic equation for the surfactant adsorption at the interface turns out to be 



dt 



7+^2^1 



In 



1- 



a 

+ r 



T 



e(^o - ipi 
T 



(3.20) 



where the last term in the brackets may be viewed as an electrostatic barrier located at 
the edge before the interface. It can be also written, by means of Eq. ( |3.9| ), as 



e(V^o - ipi 
T 



AtxI 



'0 ' 



where I = e^/eT is the Bjerrum length (about 7 A for water at room temperature). The 



kinetic equation can then be expressed, like its non-ionic parallel, Eq. (|2.14| ), in terms of 
(pQ and (f)f alone. Finally, since we have assumed negligible adsorption of counter-ions at 
the interface, we may require 

f^O, ,3.., 

The Smoluchowki equations (|3.18|) together with the Poisson equation ( |3.8| ) make a 
set of three differential equations for the three unknown functions, (j)'^{x,t), (j)~{x,t) and 
ip{x,t). Equations ( |3.9| ), (|3.19|) , ( p.20|) and ( |3.21| ) set the boundary conditions for these 
functions at the interface. If we add appropriate boundary conditions at infinity (where 
the volume fractions converge to their bulk value and the electric potential vanishes) and 
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initial conditions (say, a perfectly uniform initial state with a vanishing electric potential), 
the mathematical problem is well posed and, at least in principle, solvable. The task of 
solving this system of equations, nevertheless, seems rather formidable, though similar 



systems have recently been dealt with numerically by MacLeod and Radke [1^. Fortu- 
nately enough, in all practical cases one can avoid the elaborate mathematical treatment, 
as will be demonstrated below. 



3.4 Limiting Cases for the Adsorption 

As in Sec. 0, we are interested again in the distinction between diffusion- limited adsorp- 
tion, where the kinetics are controlled by the diffusive transport inside the solution, and 
kinetically limited adsorption, where the process is controlled by the kinetics at the in- 
terface. The kinetics inside the solution are governed now by Eqs. (|3.18| ). In order to 



identify the corresponding time scale, let us assume for simplicity that the ions have equal 
diffusivities, = D~ = D, and that the electric field is weak, etp/T <^ 1. If we add the 
two equations (|3.18|) and recall the Poisson equation (|3.8|) , the sum 0+(x,t) + (j)~{x,t) is 



found, to first order in the electric field, to undergoes a free diffusion. Since both at t = 
and t oo this sum is 0^ + (p'j^ for any x, it follows that it remains unchanged during 
the entire process. Keeping this conclusion in mind and again making use of the Poisson 
equation ( |3.8| ), we now subtract the two equations (|3.18| ) and obtain a linear diffusion 
equation for the electric potential, ip, 

dt " \ dx'^ 



^^dU-^'A- (3.22) 



This equation describes the kinetics of relaxation of the electric double layer. Its char- 
acteristic length scale is, as expected, the Debye-Hiickel screening length, and the 
resulting time scale is 

Te = {k'^D)-\ (3.23) 

More rigorously, we can calculate the asymptotic solution of Eq. ( |3.22| ), obtaining the 
following interesting time dependence, 

^(.,*^oo)==V.o(t-^)e(*-^)e- (3.24) 
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where G(t) is the step (Heaviside) function. This potential may be viewed as a "retarded 
equihbrium potential": a point x in the solution is in some sense in equilibrium with 
the interface (hence the exponential profile), but with the interface as it was some time 
{x/2kD) ago. The information propagates from the interface to the solution with the 
velocity 2k,D. From Eq. (|3.24|) it is clearly verified that the length scale is k~^, and the 
time scale is Xe as defined in Eq. (|3.23|) . 

Equation ( p.23|) states the time scale needed for the electric double layer to adjust 
and attain equilibrium with the interface. Yet, in a salt-free solution, the surfactant and 
counter-ion profiles themselves construct the electric double layer, and therefore, Te is also 
the time scale of relaxation of the surfactant profile. 

Typical values for D are about 10~^ cm^/sec, and in salt-free ionic surfactant 
solutions amounts to hundreds of Angstroms. This yields very small values for Tg (on the 
order of microseconds). The relaxation of the profile, therefore, is much faster in the case 
of ionic surfactants (without added salt) than in the case of non-ionic ones. This effect is 
due to the strong electrostatic interactions which drastically accelerate the kinetics inside 
the solution. In other words, the diffusion inside a salt-free ionic surfactant solution is an 



ambipolar diffusion ||3^ rather than a regular one. If the ion diffusivities are not assumed 
to be equal, one should replace D in Eq. (|3.23| ) with some effective diffusivity. If the 
electric field is not weak, as is practically always the case in salt-free surfactant solutions, 
then the time scale of profile relaxation will be even shorter, and the above conclusion, 
obtained for weak electric fields, will not change. 

Turning to the kinetics which take place at the interface, we treat Eq. ( |3.20| ) similar 
to Eq. ( |2.14| ) of Sec. ^ expand it close to equilibrium and find the time scale. 





= Tk^ eXp[e{lljo^eq + 'ipl,eq)/T], (3.25) 

where r^^"* denotes the kinetic time scale found in the absence of electrostatics [Eq. ( |2.17|) ]. 
As expected, the electrostatic repulsion of surfactant ions from the charged interface slows 
down the adsorption process. The sum ipo + ipi may also be written as 2ipi + {jpQ — ipi), 
where the first term expresses the slowing down due to the sub-surface concentration 
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(which is lowered because of ipi), and the second accounts for a further slowing effect due 
to the edge electrostatic barrier. If the electric field is strong, as is practically the case 
in salt-free surfactant solutions, the duration of the process may become longer by orders 
of magnitude. This can be verified experimentally when the electrostatic interactions are 
screened by added salt |T^. Using Eq. (|3.13| ) in the hmit of strong fields, together with 



Eq. (3^), we can estimate the slowing factor by 



exp[e(V'o,e, + ^i,e,)/T] exp (-^0oO • (3-26) 



This factor is typically very large. For example, in the experimental system of Ref. 

4'0,eq 



one finds na'^ ~ 10 ^, (t>teq/'Pb ~ ^^'^ l/a^ ~ 1, so the slowing factor amounts to 



about 10^ |7|. 

We see that the strong electrostatic interactions present in salt-free ionic surfactant 
solutions drastically shorten the time scale of diffusion inside the solution and drastically 
lengthen the time scale of kinetics at the interface. We expect, therefore, that ionic 
surfactants in salt-free solutions should exhibit kinetically limited adsorption. On one 
hand this conclusion greatly simplifies the mathematical treatment of the problem: we 
can safely assume that the electric double layer is in quasi-equilibrium with the changing 
interface {i.e., it obeys the Poisson-Boltzmann theory) and deal with the kinetics at the 
interface alone. On the other hand, the conclusion of kinetic control invalidates some of 
the assumptions employed by previous models. It implies that the relevant adsorption 



scenario of the two described in Sec. p.4| is the second one, illustrated in Fig. 2b. We 
recall that according to this limiting scenario, the interfacial contribution to the free 
energy, fo{4>o), retains the same shape throughout the process and reaches its minimum 
only at equilibrium. Hence, one cannot use in this case the equilibrium equation of state, 
such as Eq. ( p.l6| ), to calculate dynamic surface tensions. 

Having realized that the adsorption is limited by the kinetics at the interface, we can 
take the ion profiles and electric double layer to be in quasi-equilibrium. The Poisson- 
Boltzmann theory, therefore, can be implemented, and the kinetic equation, ( p.20| ), is 



19 



rewritten as 

+ ^ + M_2sinh-(A«)) 

(3.27) 

The mathematical problem of finding the time-dependent surface coverage, (poit), is thus 
reduced to a single integration. 



exp[(47r//a"' 



dt 



In 



1]' 



3.5 Dynamic Surface Tension in a KLA Process 

In order to relate calculated surface coverages to measured surface tensions, we need, as in 
Sec. ^, an appropriate relation between these two variables. As we have just concluded, the 
equilibrium equation of state will not do in this case. Since kinetically limited adsorption is 
considered, the contribution to the free energy from the bulk solution has the equilibrium 
dependence on the surface coverage, i.e., fluiki'Pt) of Eq. ( |3.15| ). Hence, we can write 



the free energy (or equivalently, the dynamic surface tension) as a function of the surface 
coverage alone, 

= ^ [ct>t In <Pt + (1 - <l>t) ln(l - <l^t) - </'o In <Pt - j{^{X<Pt)' + l - 1) 
+ 20+sinh-i(A0+)] - «4 - f • (3.28) 



Expression ( p.28| ) determines the behavior of the dynamic surface tension as the surfactant 



ions adsorb at the interface and the surface coverage increases. Note, that it is very 
different from the equilibrium relation ( |3.1(j| ), and hence, using the equilibrium equation 



of state to relate A7 to (pQ is invalid in this case. 

Assuming a strong electric field (A0^ ^ 1), the function A7(0^) can be shown to 
depart from the convex shape of a simple well, if P/T > 2(2 + y^) — 7.5. For such high 
values of P we should expect, therefore, an interesting time dependence of the dynamic 
surface tension, as demonstrated in Fig. 6. It should be stressed, that the curve of 
Fig. 6a is not presumed to exactly correspond to the experimental results reproduced in 
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Fig. 6b. Our claim is that an interfacial free energy of the form suggested by Eq. ( |3.28| ) 
may clearly account for the unusual, experimentally observed dynamic surface tension 
depicted in Fig. 6b. 

For values of parameters other than those chosen in Fig. 6, the interfacial free energy 
may have a non-monotonic, double-well shape. In such a case, we still expect the time 
dependence of the surface tension to schematically resemble the curve in Fig. 6b, i.e., to 
exhibit a period of almost constant tension as the system undergoes the transition from the 
first well to the second. Our current, diffusive formalism, however, cannot quantitatively 



describe the kinetics of such a process [BH, 36 



3.6 The Effect of Counter-Ion Adsorption 

An issue still to be addressed is whether such high values for the interaction constant, 
/3, are reasonable for ionic surfactants. Measurements on non-ionic surfactants yield 



typical values for P which are smaller than 4T p9|. There is no obvious reason why 
ionic surfactants should exhibit, apart from the Coulombic repulsion, also much stronger 
lateral attraction at the interface. The answer lies in one of the assumptions we have used. 
Throughout the analysis above, it was assumed that the counter-ions are surface-inactive 
and hence not present at all at the interface itself. In fact, this assumption is inaccurate 
and was taken merely for the sake of a clearer and simpler discussion. The high surface 
potentials involved in salt-free surfactant solutions should attract a certain amount of 
counter-ions to the interface. Our aim now is to calculate the correction introduced into 
the interfacial free energy as a result of the presence of a small amount of counter-ions at 
the interface. 

For simplicity, we assume the counter-ions to be in quasi-equilibrium throughout the 
system, and therefore obey the Boltzmann distribution (|3.10|) , 

Since the surface charge consists now of both surfactant ions and counter-ions, expression 
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for the surface potential should be modified, 

2T 

^0 = — sinh-i{A[0+ - (a+/a-)Vo1}- (3-29) 

In the limits of a strong field {^(pQ ^ 1) and a small density of adsorbed counter-ions 
[00 /(a~)^ ^ 00 /(a^)^], we find from Eq. ( |3.29|) that the counter- ion surface coverage 
depends quadratically on the surfactant one, 

0o-^4AV6"(0o+)2. (3.30) 

The adsorbed counter-ions introduce the following direct correction to the interfacial 
free energy as was formulated in Sec. 0, 

6{A'j)counter = t\^[T{4>o In 0o" - 0o ) " 600 ^Z-q " f^l4>o]- 

[a ) 

In quasi-equilibrium this expression reduces to 

T 

S{A'y) counter = "7 7^00 , 



SO the correction to {a'^^A'-f/T takes the form 



2 



^ , y , ,2A— <f>^{<f>ty. (3.31) 

\ / counter \ / 

In addition, the adsorbed counter-ions introduce indirect corrections into the free energy, 
through the terms e0o'0o and /^^^^^ of Eq. ( ^.2^ ). However, these two corrections turn 
out to exactly cancel each other. Note that apart from electrostatic interactions, we have 
neglected any other, short-range interactions between surfactant ions and counter-ions at 
the interface. 

Looking at expression (|3.31|) , we find that the effect of a small amount of adsorbed 
counter-ions is equivalent to an effective increase in the lateral attraction term, — /?(0o')2/2 
of Eq. (3^). The coefficient (2Aa+/a^)20{^ can be also written as / [a^Y . Taking 

a~ < a"*" ~ / leads to a correction of order 1-10, i.e., the counter-ion adsorption may 
introduce, indeed, a significant addition to the attraction term. Note, that the increase 
in (3 is sensitive to molecular details (the dimensions of the two types of ions). This fact 
may be related to the experimental observation that some ionic surfactants do exhibit the 



unusual time dependence depicted in Fig. 6b, while others do not [BS 
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4 Ionic Surfactants with Added Salt 



In many practical cases, the solution contains, in addition to the surfactant ions and their 
counter-ions, also a certain amount of dissolved salt. The effect of adding mobile ions, 
whose concentration usually exceeds that of the surfactant, is to screen the electrostatic 
interactions. Since it was found in Sec. ^ that strong electrostatic interactions drastically 
affect the nature of the adsorption process, we should expect the results in the presence 
of added salt to be significantly different from the case of salt-free ionic solutions. 

In principle, adding salt introduces two additional degrees of freedom into our formal- 
ism, namely the profiles of the salt ions and their counter-ions. This should significantly 
complicate the already elaborate problem of ionic surfactant adsorption. In order to avoid 
such complications we adopt the following simplifying assumptions: 

• The salt is surface-inactive and its concentration is much larger than that of the 
surfactant. This assumption allows us to separate the roles played by the two types 
of ions: the surfactant ions adsorb at the interface and build up the surface charge, 
whereas the salt ions form the electric double layer inside the solution. 

• The salt ions are much more mobile than the surfactant ones. Relying on this 
assumption, the kinetics of the salt ions can be neglected, and the double layer 
they form can be assumed to maintain quasi-equilibrium with the changing surface 
charge. 

• For simplicity, we assume that the salt ions and counter-ions are monovalent as well, 
or more generally, have the same valency ||39 . 



Most of the formulation given in Sec. H remains valid, in particular the equations gov- 



erning the surfactant kinetics, Eqs. (|3.18|) - (|3.20|) . The main difference is that the electric 



potential, is no longer regarded as an independent degree of freedom coupled to the 
surfactant profile, 0. Instead, it is simply given by the potential of an equilibrium double 
layer, depending on the surface coverage and salt bulk concentration alone. The high salt 
concentration we have assumed, allows us to take a screened double-layer potential in the 
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linear, Debye-Hiickel regime 



0o(i^)e" 



(4.1) 



bearing in mind that in this section, the salt bulk concentration, c<j, replaces the surfactant 
one, Cft, in the definition of k. Substituting this potential in Eqs. (|3.18| ) and (|3.19|) , the 
equations determining the kinetics inside the solution are obtained. 



'dt 
dt 



D 

D 

a 



dx \ dx 



K 



dx 



2csa? 

,2 



5oe" 



K 



dt 



(4.2) 
(4.3) 



In the last equation we have also assumed that the Debye-Hiickel screening length is 
much larger than the surfactant molecular dimension, /ta <^ 1. The equation describing 
the adsorption kinetics at the interface, Eq. ( p.20| ), remains valid as it is. 

As in the previous two sections, we are interested in the distinction between the two 
limits — diffusion-limited vs. kinetically limited adsorption. In order to find the time 
scale of diffusion, we treat the electric field as a small perturbation and seek a solution to 
Eqs. ( |4.2| ) and ( |4.3| ) which is close to the one in the non-ionic case. A procedure similar 
to the one given in Sec. ^ is now employed, yielding the following asymptotic expression. 



K.<pb<Po. 



eq 



eq 



2c,a^ 



Cb 

2c. 



eq 



2c=a^ 



3cb 
2c. 



(4.4) 



When this result is compared to its non- ionic parallel, Eq. (|2.15|) , two observations are 
to be made. The first is that Eq. ( ^.4| ) contains an additional constant term, since the 
equilibrium sub-surface concentration differs from the bulk one in the presence of an 
electric field. The second is that the diffusion time scale is slightly corrected, as expected, 
by the weak electrostatic interactions, 

2 



SO) 



1-^ 

2c, 



eq 



2c,a^ 



Scfe 
2c. 



(4.5) 



where rjf^ denotes the diffusion time scale found in the non-ionic case [Eq. (|2.16| )]. The 
correction vanishes when we take very high salt concentrations, leading to a complete 
screening of electrostatic interactions. 
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Turning to the case of KLA, we find that the expression for the time scale of kinet- 
ics at the interface derived for sah-free solutions, Eq. (|3.25|) , remains valid also in the 
presence of added salt. However, unlike the case of Sec. |^, the electric field in this case 
is weak, and the resulting slowing factor is small, so that > Tjf\ Since the time 
scales of both the diffusion inside the solution and the kinetics at the interface have been 
shown, in the presence of salt, to differ only slightly from those of the non-ionic case, we 
conclude that ionic surfactants with added salt, like non-ionic surfactants, should exhibit 
diffusion-limited adsorption. As in the case of non-ionic surfactants, this conclusion is 
well supported by experiments, as illustrated in Fig. 7. Indeed, measurements on ionic 
surfactant solutions with salt can be well fitted by theoretical curves, using the same 
schemes used for non-ionic surfactants (see Sec. 0) |12[. 

5 Concluding Remarks 

In this work we have presented an alternative approach to the problem of the kinetics of 
surfactant adsorption. One of the advantages of this approach is that the diffusion inside 
the aqueous solution and the kinetics of adsorption at the interface are not introduced as 
two separate, independent processes, but both arise from the same model. This makes the 
model more complete than previous ones, and allows us to point at the process limiting 
the kinetics of the entire system in various cases. We find the adsorption to be limited by 
the bulk diffusion in the cases of non-ionic surfactants and ionic surfactants with added 
salt, and by the kinetics at the interface in the case of salt-free ionic surfactant solutions. 
Such conclusions lead to a significant mathematical simplification of the statement of the 
problem. They are also in agreement with experimental findings. 

Another advantage is that the formulation can be readily extended to more compli- 
cated systems. We have used this to account for the kinetic adsorption of ionic surfactants 
with and without added salt. In particular, we have been able to explain the recently 
reported unusual time dependence of the surface tension in salt-free ionic surfactant solu- 
tions. In addition, our free-energy approach provides a general, straightforward method 
for calculating dynamic surface tensions from surface coverages, which does not rely on 
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equilibrium relations. This feature turns out to be essential in the case of salt-free ionic 
surfactant solutions. 

The adsorption of ionic surfactants behaves very differently in the presence or absence 
of salt: it is diffusion-limited in the former case, and kinetically limited in the latter case. 
This has been shown both experimentally and by our theory. In order to reach better un- 
derstanding of the kinetics of surfactant adsorption, additional experiments are required, 
particularly on ionic surfactants. All experiments until now have involved aqueous solu- 
tions which are either free of salt or containing high concentrations of it. It should be 
interesting to examine aqueous solutions with low salt concentrations and observe the 
crossover from one limiting behavior to another. Moreover, since the adsorption in the 
case of salt-free ionic surfactant solutions has been found to be kinetically limited, dy- 
namic surface tension measurements may be used to "probe" the actual dependence of the 
interfacial free energy, /o(0o)? on the surface coverage, 0o- This, in turn, may help explain 
the equilibrium phase behavior of ionic surfactant monolayers under various conditions 
{e.g., under compression). 

The theory presented in this work is incomplete in two main aspects. The first is that 
our diffusive formalism, as was mentioned in Sec. |^, cannot fully describe the kinetics 
in cases where an energy barrier must be overcome before the interfacial free energy 
reaches its minimum. According to the analysis of Sec. |^, actual dynamic surface tension 
measurements imply that the free energy in the case of certain ionic surfactant solutions 
probably does exhibit such a barrier. Quantitative treatment of the evolution of such 
systems should require, therefore, a more accurate (perhaps "Kramers-like" |^) theory. 
The second aspect is the lateral diffusion at the interface, whose time scale has been 
completely neglected by this work. As was mentioned in Sec. ^ cases where lateral 
diffusion is significant are encountered in practice, and a future, more complete theory 
cannot ignore its effect. 
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Figure Captions 



Fig. 1 Schematic view of the system. A sharp, fiat interface separates a dilute aqueous 
solution of non-ionic surfactants from an air or oil phase. 

Fig. 2 The two limiting cases for the time dependence of the adsorption: (a) Diffusion- 
limited adsorption (DLA) — the surface coverage, (pQ, is determined by the minimum 
of the surface free energy, /q, yet the shape of /o changes with time. The curve 
is shown at three different times, corresponding to increasing values of the sub- 
surface volume fraction, 0i: (i) at ti close to the beginning of the process, when 
01 = 5 X 10~^, (ii) at a later time t2, when (pi = 1.5 x 10~^ and (iii) at equilibrium, 
when 01 = 06 = 5 X 10~^. The energy constants are set to the (realistic) values 
a = 12T and (3 = 3T; (b) Kinetically limited adsorption (KLA) — the shape of /o 
is fixed and 0o increases with time until reaching equilibrium at the minimum of /q. 

Fig. 3 A variety of non-ionic surfactants were experimentally found to exhibit DLA. Four 
examples of dynamic surface tension measurements are given: 9.49 x 10~^M of decyl 
alcohol (open circles), as adapted from Ref. 2.32 x 10~^M of Triton X-100 



(squares), adapted from Ref. 0; 6 x lO'^M of CizEOg (triangles) and 4.35 x lO'^M 
of CioPY (solid circles), both adapted from Ref. [^. Note the asymptotic t^^/^ 
behavior, characteristic of DLA and shown by the solid fitting lines. 

Fig. 4 (a) The dependence between the surface tension and surface coverage in a DLA case 
(the parameters are the same as in Fig. 2); (b) The schematic time dependence of 
the surface tension expected from a dependence A7(0o) such as in (a): an initial 
slow change, then a rapid drop and finally a relaxation towards equilibrium. 

Fig. 5 Schematic view of the ionic salt-free system. A dilute solution containing surfactant 
ions and their counter-ions has a sharp, fiat interface with an air or oil phase. 

Fig. 6 (a) The dependence between surface tension and surface coverage in the KLA case, 
according to Eq. ( p.28|) , where the values taken for the parameters are: a'^ = 12 A, 



3.6 X 10 ^, a = 11. 5T and (3 = 9T; (b) The experimentally observed dynamic 
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surface tension of a salt-free 3.5 x lO^'^M SDS solution against dodecane, adapted 
from Ref. [|l^. The time dependence is schematically consistent with a non-convex 
dependence A7(0q ) such as the one in (a). 

Fig. 7 Ionic surfactants in the presence of added salt were experimentally found to exhibit 
DLA. Three examples are given: Dynamic surface tension of 4.86 x 10~^M SDS 
with O.IM NaCl against dodecane (open circles and left ordinate), adapted from 
Ref. [mi; Dynamic surface tension of 2.0 x IQ-'^M SDS with 0.5M NaCl against air 
(squares and left ordinate), adapted from Ref. Surface coverage deduced from 



Second Harmonic Generation measurements on a saturated solution of SDNS with 
2% NaCl against air (solid circles and right ordinate), adapted from Ref. [JT|. Note 



the asymptotic t behavior, characteristic of DLA and shown by the solid fitting 
lines. 
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